function f = fstif(u,t)
f = -1e6 * (u - cos(t)) - sin(t);
end

